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Abstract 

J-PET analysis framework is a flexible, lightweight, ROOT-based soft¬ 
ware package which provides the tools to develop reconstruction and cal¬ 
ibration procedures for PET tomography. In this article we present the 
implementation of the full data-processing chain in the J-PET framework 
which is used for the data analysis of the J-PET tomography scanner. 
The Framework incorporates automated handling of PET setup parame¬ 
ters’ database as well as high level tools for building data reconstruction 
procedures. Each of these components is briefly discussed. 


1 Introduction 

Positron Emission Tomography (PET) is one of the most accurate, non-invasive 
techniques used in modern medical diagnostics to produce the metabolic image 
of patient by reconstructing the spatial distribution of a radiopharmaceutical 
administered into the patient’s blood. The signals coming from the patient’s 
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body in the form of gamma quanta are recorded by scintillator-based tomogra¬ 
phy scanners. 

The Jagiellonian PET (J-PET) collaboration develops a novel concept of 
a Time-of-Flight PET (TOF-PET) apparatus based on plastic scintillators, 
which largely improve the field of view and the timing resolution of the TOF 
method urn- In order to achieve this goal, the J-PET detector consists of 
large strips of plastic scintillators read out at both ends by photomultipliers [6]. 
The photomultiplier signals are sampled and recorded by dedicated Front-End 
Electronics (FEE) specifically designed for this project T[[8]. 


2 Technology 

J-PET analysis framework has been developed mainly in the C-|—b program¬ 
ming language applying the object-oriented approach. Additionally, small parts 
of the code are written in Python and Perl. The C++11 standard and Boost 
Libraries have been used. Boost Unit Test Framework [5] has been chosen to 
construct a set of automatic test to assure the quality of the code. The package 
is based on the Open Source library ROOT m which is a de facto standard in 
most of the modern high-energy physics and hadron physics experiments. The 
documentation of the code is automatically generated using dedicated Doxy- 
gen program m- Besides, several external libraries are used: tinyxml m and 
libconfig++ H3]. The database handler module responsible for the communica¬ 
tion with the PostgreSQL database uses PostgreSQL Application Programming 
Interface (API) T4(. The SQL queries have been enclosed in the PostgreSQL 
functions which are remotely called from the framework code. 


3 Framework architecture 

The basic architecture of the J-PET framework was already explained in ref. Ca¬ 
in this section we shortly review the main idea and explain in more detail the 
integration of the framework with the parametric database. 

The data reconstruction in modern PET tomography is a complex multi¬ 
stage process including several low and high level reconstruction algorithms and 
calibration procedures, which needs a special computing support m- Also, 
some monitoring procedures are required to be implemented. The basic concept 
of the J-PET Framework is the decomposition of analysis chains into series of 
standardized modular blocks. Each module corresponds to a particular comput¬ 
ing task, e.g. reconstruction algorithm or calibration procedure, with a defined 
input and output methods. The processing chain is built by registration of cho¬ 
sen modules in the JPetManager, which is responsible for the synchronization 
of the data flow between the modules. This approach has several advantage e.g. 
it permits to quickly interchange modules and to create a processing chain that 
can be used for different experimental setups, e.g. in the J-PET project, the 
substitution of the lowest level processing modules allows to analyse both the 
raw data delivered by the full Data Acquisition system (DAQ) but also data 
registered in the diagnostic setup using a digital oscilloscope (see Fig. 0. 

In addition, the framework is well suited to develop procedures in paral¬ 
lel, by taking advantage of the standardize module structure and predefined 
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Figure 1: Scheme of an example data processing chain with the J-PET Frame¬ 
work. 


input/output. This feature is particularly important in the prototype phase 
of the project, where the several reconstruction procedures are developed and 
tested at the same time by different subgroups. 

3.1 Data processing 

The J-PET Analysis framework aims at providing the possibility to implement 
reconstruction procedures to users even with little programming experience 
through a convenient and simple API abstracting the technical details. Cre¬ 
ation of a simple analysis stage requires providing a module class which basi¬ 
cally comprises three methods only: one to be called prior to data processing, 
one after processing, and a method called for each event in the dataset. Any 
non-standard processing can be achieved, however, simply by overwriting more 
specialized methods inherited from the base module. The elementary pieces of 
recorded data are enclosed in objects containing not only the experimental and 
reconstructed information such as recorded time signals from Time-to-Digital 
Converters (TDC), but also pointers to objects representing parameters of all 
the elements of the DAQ system related to each signal. These pointers are im¬ 
plemented using the TRef mechanism from the ROOT framework which allows 
them to be persistently stored and retrieved from ROOT files. Therefore, the 
set of parameters from the database is downloaded only once before the process¬ 
ing and then stored along with data files which allows to access each necessary 
parameters related to a certain signal with 0(1) operations at any time. For 
every class a list of handful methods defines an interface that makes the access 
to data easy. The relation between the parametric objects reflexes the logical 
or physical connection between the objects they represent, e.g. JPetPM object 


3 











describing the proprieties and setting of a photomultiplier, points to JPetScint 
object, which corresponds to a scintillator module to which this photomulti¬ 
plier is connected. The information about all the relations is encoded in the 
parametric database and then is restored while generating parameter objects. 

A typical processing scheme with the Framework is split between several 
modules, each analysing a set of data objects and combining them into more 
complex structures as shown in Fig. |T] The analysis starts with a set of raw 
TDC signals and uses them to reconstruct shapes of electric signals from pho¬ 
tomultipliers. Next, pairs of signals corresponding to a single photon hit in a 
scintillator strip are identified and used to build hit objects. Finally, pairs of 
hits from a single annihilation event are merged into Event structures which in 
turn constitute the input to image reconstruction procedures designed for the 
J-PET (TTlITTil . Each of these steps is performed by a separate module which 
allows for easy exchange of the procedures in case of change in the DAQ system. 
For instance, tests described in [6! were performed with a Serial Data Analyzer 
providing dense signal sampling with the dedicated DAQ, but the respective 
data can be processed by the same framework with an exchange of the lower- 
level modules while the higher-level processing modules are common for all DAQ 
chains (see Fig. |T|). 

3.2 Connection with the parametric database 

During the PET examination all information about detector setup and experi¬ 
mental conditions must be stored. In the J-PET project the dedicated paramet¬ 
ric database is used for this purpose [20]. It uses PostgreSQL object-relational 
database management system. The database contains records about hardware 
settings, calibration constants, initial setup parameters, alignment of the detec¬ 
tors and software configurations. For each single measurement all these settings 
must be saved and available at any time for analysis and data processing. The 
communication between the parametric database and the analysis framework is 
done by a dedicated module (DBHandler) [21]. It works as an interface handling 
all connection operations and hiding the details of the database structure from 
users. The Parameter Manager generates necessary information for analysis and 
reconstruction modules, based on the data provided by the DBHandler, in the 
form of collections of parametric objects describing geometry, scintillators, pho¬ 
tomultipliers, voltage etc. The Parameter Manager is also responsible for saving 
and reading the subsets of parameter objects in ASCII or ROOT file formats. 
The information retrieved from the database is stamped with the run number, 
which uniquely determines a given measurement condition set. The whole pro¬ 
cess is fully automated. The user has to give a valid run number identifier and 
the appropriate set of parameters will be automatically created and loaded in 
its analysis module. 


4 Summary 

In this article we have presented a software framework for the PET tomography 
data analysis and reconstruction, that has been developed as a part of the J-PET 
project. The package is ROOT-based and it is inspired by several solutions from 
the hadronic and high energy physics experiments |22U'25j . The JPet Framework 
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provides a set of programming tools and predefined data objects that helps in 
the implementation process even for unexperienced users. As it was shown in 
the examples, the proposed modular architecture permits to efficiently handle 
different data sources, DAQ setups and detector geometries by exchanging the 
low-level modules in the processing chain. Finally, the J-PET Framework can 
provide information on the measurement conditions, geometries and hardware 
setups as a set of parameter objects generated based on the database registers in 
a fully automated way. The first full analysis chain starting from the ’’raw” data 
delivered by scanner data acquisition system up to the reconstructed set of Line 
of Response has been implemented. Currently, the full signal reconstruction in 
scintillator methods based on the Compressive Sampling |26j and the library of 
signals m are being incorporated. 
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